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Abstract 

Over the last 15 years, the supernova community has endeavoured to directly identify progenitor stars of 
core-collapse supernovae discovered in nearby galaxies. These precursors are often visible as resolved stars 
in high resolution images from space and ground based telescopes. The discovery rate of progenitor stars 
is limited by the local supernova rate and the availability and depth of archive images of galaxies, with 
18 detections of precursor objects and 27 upper limits. This review compiles these results (from 1999 - 
2013) in a distance limited sample and discusses the implications of the findings. The vast majority of the 
detections of progenitor stars are of type II-P, II-L or lib with one type lb progenitor system detected and 
many more upper limits for progenitors of Ibc supernovae (14 in all). The data for these 45 supernovae 
progenitors illustrate a remarkable deficit of high luminosity stars above an apparent limit of logL/L0~ 
5.1 dex. For a typical Salpeter initial mass function, one would expect to have found 13 high luminosity and 
high mass progenitors by now. There is, possibly, only one object in this time and volume limited sample 
that is unambiguously high mass (the progenitor of SN2009ip) although the nature of that supernovae is 
still debated. The possible biases due to the influence of circumstellar dust, the luminosity analysis, and 
sample selection methods are reviewed. It does not appear likely that these can explain the missing high 
mass progenitor stars. This review concludes that the community’s work to date shows that the observed 
populations of supernovae in the local Universe are not, on the whole, produced by high mass (M )> 18Mq) 
stars. Theoretical explosions of model stars also predict that black hole formation and failed supernovae 
tend to occur above an initial mass of M ~ 18Mq. The models also suggest there is no simple single mass 
division for neutron star or black-hole formation and that there are islands of explodability for stars in the 
8 — 120Mq range. The observational constraints are quite consistent with the bulk of stars above M > 18Mq 
collapsing to form black holes with no visible supernovae. 

Keywords: (stars): supernovae: general - stars: Wolf-Rayet - stars: massive - stars:evolution - (stars): 
supergiants 


1 Introduction 

The link between massive progenitor star and the type 
of core-collapse event it produces at the end of its nu¬ 
clear burning live is the fundamental piece of informa¬ 
tion that underpins our understanding of the physi¬ 


cal processes involved in stellar explosions, (see Langer 


2012L for a recent theoretical overview). Since the end 


of the 1990’s the supernova community has been exten¬ 
sively searching publicly available archives for high res¬ 
olution and deep images of nearby galaxies which host 
supernovae (SN). The existence of these pre-explosion 
images has allowed direct identification of the progen- 


novae (e.g. Van Dvk et al. 

2003b 

: Smartt et al. 2004: 

Maund et al. 20oll Li et al 

2006 

). Many of these, as 


individual events, are quite compatible with stellar evo¬ 
lutionary model predictions. The fact that numerous 
red supergiants are now observed as the stellar end 
points that produce the most common class of type II 
SN is a validation of modern stellar structure theory. In 
2009, the number of detections of progenitor stars had 
reached a po i nt tha t analysis of samples was possible 
( Smartt et al. 20091) and a review of what these con¬ 
straints me ant for the_Jield in the broader sense was 
warranted ( Smarttl 120091) . Since then there has been 
further very significant progress and results which are 
more consequential than just an increase in sample size. 
This review will consider work done in the field since 
2009, and review the results as a whole in the context 
of the SN population we observe in the Local Universe 
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as quant if ied by , for example, Leaman et al.1 ( 2011 ). and 
Li et al. 02011 ). Some of the highlights in the last 5 
years are the significant increase in the sample size, nine 
further detections of progenitors, evidence for the disap¬ 
pearance of progenitor stars many years after explosion 
and detections of possible binary companions. There 
has also been active discussion in the literature concern¬ 
ing the analysis methods, particularly in the treatment 
of extinction toward progenitors and how that may af¬ 
fect the luminosity and mass estimates. There are even 
some progenitor stars with multi-epoch data that allow 
pre-explosion variability to be probed, and as data sets 
grow then lightcurve monitoring of progenitors may be¬ 
come more common place (e.g. Elias-Rosa et al. 20091: 
Szczvgiel et al • 2012 : Fraser et all[20l3 h 
In an att e mpt to set a well defined sample , 


Smartt et al.l ( 20091) (and later lEldridge et al.1 12013 ) 
defined criteria for selecting a time and volume limited 
sample of SNe for which searches for progenitors 
were feasible : all CCSNe in galaxies with recessional 
velocities V v - lr < 2000kms _1 (corrected for Virgo in¬ 
fall, which corresponds to a distance of d < 28Mpc 
for Hq = 72kms -1 Mpc -1 ). This distance limit was 
imposed from practical experience of analysing HST 
archival image data and the difficulty in retrieving 
photometry of resolved individual massive stars beyond 
this distance. Additionally, as galaxy number count 
increases as the distance cubed, the relative number 
of galaxies with high quality archival imaging rapidly 
decreases. Hence the probability of finding coincidences 
of SN discoveries combined with existing pre-explosion 
images of the galaxies also falls off quickly. As we will 
see later in this review the SN which yield the most 
restrictive information on their progenitor stars have 
tended to be significantly closer than 28 Mpc. Two 
groups have worked extensively in this area, a Cali¬ 
fornia based t eam (who started the ir concerted search 
with HST in IVan Dvk et al. 19991) , and a UK based 
(initiated with their first paper: 


Smartt et al 


but others have also made significant contribu¬ 


tions (e. g. iGal- Yam et al' 2007 : Kochanek et al. 2012 : 


Prieto et all 20081 ). The IJK based group have focused 


their attention on this volume and time limited survey 
to get high-resolution, deep imaging of all SN that have 
useful archival imaging so that the precise SN positions 
can be astrometrically placed on the pre-explosion im¬ 
ages. The Van Dyk led team have also pursued a similar 
strategy and bot h teams have used combinations of 


2011: Van Dvk et al. 2014) g 

round-based AO imaging 

(e.g. Fraser et al. 2011, 2012| 

Van Dvk et al. 2012a) or 


good image quality natu ral seeing images (e.g. iLi et al 


20051 : Fraser et ah 20141) to locate the SN positions to 
accuracies of ~30-50 milliarcsec on the pre-explosion 


frames. 


2 Volume and time limited survey data 

The definition of a time and volume limited survey 
has significant advantages in determining the overall 
properties of supernova progenitor systems as a func¬ 
tion of initial mass. There have been some remark¬ 
able discoveries of individual stars before explosion, 
with tight ph ysical constraints on their luminosity and 
temperature ( Smartt et al.l2004 ; lVan Dvk et all2012bl ) 
some of which have been shown to have d isappeared 
( Van Dvk et al. 2013 : Maund et al. 2014bl) . However 
the most interesting aspect of this area is now the sta¬ 
tistical sample of more than a dozen unambiguous di¬ 
rect detections and many restrictive limits to produce 
a defined survey sample of 44 objects. Objects outside 
this distance limit have also contributed to our under¬ 
standing, such as t he massive progenitor of SN2005gl 
at 60Mpc ( Gal-Yam et al. 2007 ) but the relative fre¬ 
quency of such an event needs to be put in context with 
the fixed volume sample. 

This review extends the time and volume lim- 

(|2009h 


ited sample introdu ced by iSmartt et al.1 ([2009) and 
Eldridge et al. ([20131 ) for the type II and Ibc progen¬ 
itors to the end of 2013. The discussions is extended to 
include all other progenitors within this survey sample 
including the lib and Iln types and some beyond. The 
closest SNe provide the best opportunity for studying 
their progenitors, carrying out multi-wavelength mon¬ 
itoring until very late times. It is likely that we are 
not missing a large fraction of SNe in galaxies closer 
than ^12 Mpc ( Botticella et al. 2012 ) as the starforma- 
tion rates of these galaxies are in quite good agreement 
with the measured SN rate for a a lower mass limit of 
8M 0 for core-collapse (although iHoriuchi et al.1 l201ll 
suggests that at higher redsliift there is a discrepancy). 
These closest SNe offer the best opportunities for pro¬ 
genitor detection, but as we can’t influence the SN rate 
nor the discovery rate (assuming we are 70-80% com¬ 
plete) the only way to increase the progenitor discovery 
rate is through patience and time. The following sec¬ 
tions discuss the results published to date for the dif¬ 
ferent subtypes of progenitors and Section[3] presents a 
discussion of these empirical results in the context of 
massive stars, the initial mass function, and the local 
SN population. 

2.1 Type II supernovae 

As illustrated in Smartt et al.l ( 20091) . the chance of 
having an image of a nearby galaxy with one of the 
three main imaging cameras onboard HST (WFPC2, 
ACS or WF3) for any nearby SNe (within T4ir < 
2000kms -1 ) is about 25%. If multi-colour imaging is 
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Figure 1. Analysis of the progeni tor of SN 2 012aw . Upper: The 
SED fit of a DUSTY model from IKochanek et alJ l| 2012f) to the 
detect i ons in the optical and N IR filters presented in lFraser et al.l 
J~20 1 2f ): [~Van Dvk et al.l ll2012al ) and the mid-IR Spitzer limits. The 
black curves show the SED model which is made up of the three 
components of scattered photons, direct photons and dust emis¬ 
sion radiation (the red dashed curves from left to right respec¬ 
tively). The blue line is the unobscured SED of the star itself. 
This is for a dust temperature of T ( j = 1000 K and a graphitic 
dust composition (see the original reference for a similar sili¬ 
cate fit). Lower: Comparison of the luminosity estimates of the 
three s tudies. The point with error bar is from |Van Dvk et alJ 
(2012a ), dashed contours fro m iFraser et al.1 ll2012ll and red locus 
from IKochanek et al.1 ll2012h . illustrating the importance of con¬ 
sidering the treatment of both dust absorption and emission. As 
discussed in Sectioi j2.1.2l the proper inclusion of dust treatment 
lowers the luminosity estimates substantially. Reproduced from 
Figures 5 and 6 in “On Absorption by Circumstellar Dust, with 
the Progenitor of SN 2012aw as a Case Study by Kochanek et al. 
2012, ApJ, 759, 20. 


available then these colours can be fit with either ob¬ 
served supergiant colours or model atmospheres to de¬ 
termine approximate spectral types or effective tem¬ 
peratures, bolometric corrections and hence bolomet- 
ric luminosity. An alternative, but equivalent approach 
is to use stellar evolutionary models that have as¬ 


sociated model atmosphere spectral energy distribu¬ 
tions (SEDs) for the surface characteristic T e g, and 
logL (or log g) and synthetic photometry calculated 
from these spectra. There have been variations on 
these two methods from the early days of the first 
red supergiant detection of the progenitor of SN2003dg 
( Van Dvk et al. 2003b : Smartt, et al.lI20 q 3 ) through to 
the latest which have extended spectral colours in 


the near infra-red dFraser et al. 2012 : Kochanek et al 


2012! : Van D vk et al. l2012af> . Detections of red super¬ 


giants in JHK -bands are quite powerful indicators of 
the bolometric flux, since extinction uncertainties are 
lower and the bolometric corrections are less uncer¬ 
tain ( Levesque et al. 20061 ). Deep enough near infra-red 
pre-explosions images of nearby SN are still the excep¬ 
tion rather than common place (although limits can 
set restrictions on su per-AGB progenitors for example; 


Eldridge et al, 2007) 


Early results indicated that the detections of the pro¬ 
genitors of II-P SNe were red supergiants, with esti¬ 
mated luminosities in the range T3< log L/L 0 <5.0 
(see review of results in Smartt! l2009h . which would 
mean zero age main-sequence (ZAMS) masses in the 
region 7-15 Mq. Since then, improvement of the anal¬ 
ysis methods have provided updated estimates of the 
stellar parameters both for previously published detec¬ 
tions and new discoveries. Efforts have been made to 
systematically include the MARCs stellar models and 
consistently apply synthetic photometry to reduce com¬ 
parative errors in colour corrections and hom ogenise t he 


analysis (for more i n dep th discussion see iFraser et al 


20111 : Dayies et al. 120131 ). Observing the SN progeni¬ 


tor position 3-4 years after explosion to check that the 
identified stellar progenitor has disappeared now allows 
difference imaging techniques to be applied in order 
to i mprove the pre c ision o f the pre-explosion p hotom- 
etry. Maund et al.1 ( 2014bl) and Maund et al.l ( 2014al) 
discuss this in detail and the case for the disappear- 
ance of the progenit or of SN2011dh is presented in 


Van Dvk et al. (I2013T) . The treatment of the extinc¬ 


tion to ward the progeni t ors ha s progressed and the pa¬ 
per of Kochanek et al. ( 20121) highlights the need to 
consider interstellar medium (ISM) and circumstellar 
medium (CSM) C SM extinction separately and consis¬ 
tently. ICIroh et all ( 2013 d) have considered rotating and 
non-rotating models showing that initial rotation of the 
progenitor introduces a further uncertainty on the pro¬ 
genitor mass estimate as the core luminosity changes. 
They suggest that any particular mass may be uncer¬ 
tain to within 4-5 Mq but the overall mass range derived 
for type II-P pr ogen it ors is q uite similar to that orig¬ 
inally derived in Smartt et al. ( 20091) . All of this work 
has put the estimates of stellar parameters on a firmer 
footing with some changes in the stellar T e g and logL 
for various progenitor stars. However it is fair to say 
that the improved estimates of stellar luminosity (and 
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hence mass) are not significantly different to the orig¬ 
inally publish ed work. The mai n result from the data 
up to 2009 ( Smartt et al.ll2000l ) still holds : type II-P 
supernovae come from red supergiants for which the 
lowest mass estimated is around 8M 0 and the highest 
mass progenitor has a luminosity of L ~ 5.OL 0 (corre¬ 
sponding to a ZAMS of 15-18M 0 ). 

The most recent results, and in the opinion of the 
author, the most reliable estimates of the the physical 
parameters of detected progenitors of type II SNe are 
listed in Table[T]with the references to the relevant pa¬ 
pers. Where new and updated estimates are used, refer¬ 
ences to the original detections are provided with com¬ 
ments. This review will not discuss all the detections 
individually, and the differences between the measure¬ 
ment and analysis methods. However two SNe are taken 
here as examples for a more in depth discussion as case 
studies. The issues either in updated stellar luminosi¬ 
ties or the disagreement between analysis methods are 
discussed for SN2004A and SN2012aw. 

The two SNe which fell on compact, and likely co¬ 
eval clusters (SN2004dj and SN2004am) in NGC2403 
and M82 are not considered here, although the turn¬ 
off mass estimated for both is in the regime of the 
moderate to low masses that are derived for II-P pro¬ 
genitors. Hence their inclusion would support the re- 
sul ts discussed below. S N1999ev is also not included 
Maund et al.l ( 2014bl) sho wed that the pro g enitor 


as 


object identified originally bv lMaund fc Smartt ( 2005) 
was likely a stellar cluster . Fina lly SN2003ie, which was 
included in lSmartt et al.l ( 20091) is no longer considered 
as with the higher extinction estimates used later in this 
analysis, the luminosity limit is not useful nor restric¬ 
tive in any way. 


2.1.1 SN2004A 

This SN progenitor is taken as an example were the 
original estimate of luminosity has a significant dis¬ 
crepancy with a new measurement. As the later mea¬ 
surement is based on new data after explosion, with 
difference imaging applied it is worth discussing as 
a an exampl e how sign ificantly resul t s cou ld change. 
Smartt et al. ( 2009h and lHendrv et al. ( 20061) originally 
estimate d a logT/L^=4.5±0 .3. With deep follow-up 
imaging Maund et al.l ( 2014bi) showed the progenitor 


object identified has disappeared and suggested a sig¬ 
nificantly higher value of log L/L 0 =4.9 ± 0.3. The es¬ 
timate of the magnitude of the progenitor in t he HST 
F814W pre -explosion image is sim ilar in iHendrv et al 


( 20061) and Maund et al. ( 2014bl) . illustrating the dif¬ 
ference in final quoted luminosity arises in the analysis 
methods. The methods differ in that a high er extinc¬ 
tion of E(B — V) = 0.16 was employed bv lMaund et al.1 
( 2014bl ). together with application of the MA RCS 
model atmospheres. In his PhD thesis work, Fraser 
( 201ll) undertook a re-analysis of the type II-P pro¬ 



Wavelength (A) 



Wavelength (A) 


Figure 2. The SED fit of the progenitor of the yellow super¬ 
giant progenito r of SN 2011dh from two independent an alyses of 
I Van Dvk et al.l (1201 lh (upper) and iMaund et ahl (1201 if) (lower). 
The pre-explosion HST data points (solid black with errors) are 
impressively fit with a single stellar source and an ATLAS model 
atmosphere for T g r = 6000 K and log g = 1.0 in both analyses. 
Reproduced from Figure 3 in “The Progenitor of Supernova 
2011dh/PTFlleon in Messier 51 by Van Dyk et al. 2011,ApJ, 
741, L28 and Figure 2 in “The Yellow Supergiant Progenitor of 
the Type II Supernova 2011dh in M51 by Maund et al. 2011, 
ApJ, 739, L37 . 


genitor sample. As an experiment, he suggested that 
adding in extinction values that are seen toward the 
red supergiant population of local group galaxies M31, 
the Large Magellanic Cloud (LMC) and the Small Mag¬ 
ellanic Cloud (SMC) would be appropriate to test if the 
typical extinctions estimated for SN progenitors are sys¬ 
tematically lower th an those observ ed toward existing 
red supergiants. lSmartt et al. ( 20091) had also discussed 
this method for adding extra ad hoc extinction to those 
progenitors for which foregrou n d esti mates only where 
available. While ISmartt et al. (2Q0§) had adopted an 
additional Ay ~ 0.3. iFraseil (1201ll) adopted an addi¬ 
tional Ay ~ 0.53. Hence there a re three estimates of 
1 okL/Lm= 4. 5, 4.8 and 4.9 f rom Smartt et al. ( 20091) : 
Fraser ( 2011 ): Maund et al. ( 201-llJ ). each with an un¬ 


certainty of ±0.3 dex. This is probably the most extreme 
case of a discrepant result due to different methods, al- 
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Table 1 List of detections of type II progenitors. The comments in the Reference column give guidance as to the source of the stellar 
progenitor data for the luminosities (logL/L©). The initial masses of model stars that end with these luminosities are given in the last 
two columns : M( S,G) are the masses (in Mq) from the STARS and Geneva rotating models (which have similar final luminosities); 
M( K) are from the KEPLER models (see Section fT 1 1 for details). The values for the initial masses from KEPLER for 2003gd, 2005cs 
and 2009md are uncertain since there are no KEPLER models evolved to end points at this mass and luminosity. They are based on 
extrapolation, and quoted in parentheses. The errors for M(K) can be assumed to be the same as those quoted for M(S,G). 


SN 

Type 

log T eff 

err 

log L 

err 

Reference 

M(S,G) 

M{ K) 

SN2003gd 

II-P 

3.54 

0.02 

4.3 

0.20 

From 

(1), similar to (2) 

7-i 

(8) 

SN2005cs 

II-P 

3.55 

0.05 

4.4 

0.20 

From 

(3), similar to (5) and (6) 

8 1\ 

(9) 

SN2009md 

II-P 

3.55 

0.01 

4.5 

0.20 

From 

(11) 

9±i 

(10) 

SN2006my 

II-P 

3.55 

0.10 

4.7 

0.20 

From 

(3), similar to (4) and (7) 

io ti 

12 

SN2012A 

II-P 

3.58 

0.05 

4.7 

0.10 

From 

(13) 

ioi2 

12 

SN2013ej 

II-P 

3.57 

0.04 

4.7 

0.20 

From 

(16) 

10±2 

12 

SN2004et 

II-P 

3.56 

0.05 

4.8 

0.20 

From 

(4) and (17) 

12^3 

13 

SN2008bk 

II-P 

3.64 

0.10 

4.8 

0.20 

From 

(3), similar logL/L 0 =4.6 in (8) 

12^3 

13 

SN2004A 

II-P 

3.59 

0.04 

4.9 

0.30 

From 

(3). see Sect|2.1.1l 

13±| 

14 

SN2012aw 

II-P 

3.56 

0.04 

4.9 

0.10 

From 

(14). see Sect|2.1.2| 

13^2 

14 

SN2009hd 

II-L 

3.72 

0.15 

5.0 

0.20 

From 

(12) 

15±1 

16 

SN2009kr 

II-L 

3.68 

0.03 

5.1 

0.25 

From 

(9), similar results in (10) 

16±| 

18 

SN2012ec 

II-P 

3.53 

0.03 

5.1 

0.20 

From 

(15) 

16±s 

18 


(I) 

(5' 

Maund & Smartt (2009): (2)Smartt 

et al. (2004); (3lMaund et al. (2014b); (4lCrockett et al. (2011); 

Smartt et al. (20091): (6) 

Maund et al. 

(2004: (7]Li et al. (2007); (SlMattila et al.l (20081): (9)Fraser et al. 

(2010); (10lElias-Rosa et al. 

(20li" (11 

Fraser et al. (2011): (12)Elias-Rosa et al. (2011) : note the T e ff and 

log 

fro 

L/L 0 values are upper limits. (13]Tomasella et al. (2013) : note the T e tr was an assumption, and not derived 

m multiple colours: (T4) Kochanek et al 

(20121): (15lMaund et al. (120131): ( 16)Fraser et al. (2014) : note the XI// 


was an assumption, and not derived from multiple colours; (17) Fraser et al. (in prep) 


though formally, the uncertainties overlap. The value of 
logL/L 0 =4.9 ± 0.3 is reported in TablcjU and it serves 
to illustrate the point made earlier that results from 
improved analysis methods should be adopted for indi¬ 
vidual objects were possible. However it will be shown 
that the luminosity shifts do not affect the important 
overall trend of a lack of high mass stars. 


2.1.2 SN2012aw 

The second specific case study example is the progeni¬ 
tor of SN2012aw. This was a very well studied type II- 
P SN in M95 (d = lOMpc) discovered within 1-2 days 
of explosion and followed well into th e nebular phase 


(Bose et al.l2013t Dall’Ora et al 


(12012al ) and IFraser et all (2012) both identified the 


2014). IVan Dvk et all 


same progenitor object as the culprit star with detec¬ 
tions in HST filters F5551T, F8141V and in near infra¬ 
red ground based images through filte rs J a nd K. The 
galax y was also imaged by the SINGS ( Kennicutt, et al.1 
2003) Spitzer survey but no detection of a single 


point source w a s pos sible, a s discussed in deta il in 
Kochane^e^aL (2012). Both Fraser et al. ( 20121 ) and 


Van Dvk et al.1 ( 2012al) carried out a similar analy¬ 


sis and suggested the progenitor star was quite lumi¬ 
nous and relatively high mass in the range 14-26M 0 
and 15-2OM 0 respectively. They also suggested that 
the progenitor suffered significant reddening. Since the 


SN itself suffered low line of sight extinction, this was 
taken as evidence that the progenitor reddening was 
due to CSM dust which was destroyed in the explo¬ 
sion. Up to this point in the history of progenitor 
analysis, all treatments of reddening towards progen¬ 
itors had assumed that applying interstellar laws and 
relations were applicable also to circumstellar dust. 
However iKochanek et al.l ( 2012 ) presented an alterna¬ 
tive analysis in which they discussed the three con¬ 
sequences of such assumptions. Attempting to model 
CSM dust with ISM laws ignores dust emission in the 
near-IR, assumes (wrongly) that all scattered photons 
are lost to the observer and do not treat the (probable) 
dust composition corr ectly. In the case of SN2012aw , 
Kochanek et al. ( 20121) argu ed that lFraser et al.1 ( 20121) 


and Van Dvk et al. ~ i|2012al) had significantly overes¬ 
timated the lu minosity of the prog enitor. Using the 
DUSTY code of lvezic fe Elitzur ( 19971) to compute how 
photons are scattered, absorbed and reemitted by a 


dusty region enclosing a stellar source. IKochanek et al 


( 2013 ) calculated a model to fit the detected magni¬ 
tudes and the upper limits from the mid-IR (see FigQ]). 
They suggested a significantly lower luminosity progen¬ 
itor, a stellar temperature of T e g ~ 3600^2oo K and an 
upper mass limit of Mzams < 15M 0 . The amount of 
CSM gas, due to the progenitor’s stellar wind was also 
constrained by the x-ray and radio detections. The mass 
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loss estimate of M 10~ 5 M 0 yr _1 is was roughly con¬ 
sistent with a star of this temperature and luminosity. 

This progenitor is important as the discussion over 
the CSM reddening has implications for estimates of 
the most massive stars that can explode as bright su¬ 
pernovae, as we shall see in Section!?] 


2.1.3 Upper luminosity limits 

While TablcQ] lists the thirteen type II SN that have se¬ 
cure detections of their progenitors in the period from 
1999-2013 (inclusive), there are a further 13 type II SNe 
which have pre-explosion HST data (or alternative high 
quality ground based images) in which no progenitor is 
detected. In these cases the progenitor falls below the 
detection limit of the images and no information on the 
colour and effective temperatures of the stars can be 
gained a priori. Any individual event does not there¬ 
fore add a lot of new information or insights compared 
to the large number of detections now available. How¬ 
ever the statistical significance of the number of non¬ 
detections does provide a constraint on the mass and 
luminosit y range of the proge nitor po pulation as dis¬ 
cussed in Smartt et al.1 ( 20091) and in Kochanek et al 


( 20081) . Since the Smartt et al. (2009) study, there 
have been 3 more non-detections that can boost the 
statistical significance of the sample further. These 


are SN2009H (Li et ali 2009), SN2009N (iTakats et al 


20141) and SN20 0 9ib (ITakats et al.ll2015l) . As discussed 
in Smartt et al] (2009 ) and previous papers s uch as 
Smartt et al. ( 20031) and Van Dvk et al.l ( 2003a ). if one 


assumes that the progenitors of these type II (mostly 
II-P) SNe are red supergiants then a luminosity limit 
can be calculated for this effective temperature regime. 
An /-band like filter sampling the wavelength region 
between 7000-9000A is particularly powerful since the 
stellar spectral energy distribution peaks in this region 
and the bolometric correction change means the lumi¬ 
nosity limit is roughly constant across the K to M-type 
supergiant colours and temperatures. The major uncer¬ 
tainty is the extinction to adopt toward the progenitor 
and different methods have been proposed. 

Initial work in this area simply adopted the estimated 
line of sight extinction toward the SN as appropriate to 
apply to the progenitor. The extinction could be esti¬ 
mate d from the spectral slope of the S N. (iSmartt et al 


2011 


2003) the absorption of Nai lines (IPoznanski et al 


120121) or from three colo ur photometry of th e 


surrounding stellar population (iMaund et al.1 2014b). 
There are problems with all of these estimates, but even 
more important is that this line of sight extinction may 
either not be directly applicable to the progenitor (as in 
the case of the surrounding stellar population) or may 
be simply a lower limit since it does not take into ac¬ 
count circumstellar medium extinction in the locality of 
the progenitor star (see the discussion in Section !?. 1.21) . 


It is quite possible, and indeed probable that CSM dust 
(within a radius of approximately 10-100/?*, where /?* is 
a typical red supergiant (RSG) radius of a 500-1000R©) 
is destroyed during the UV flash of the SN. The discov¬ 
ery of a dust enshrouded progenitor, and a re latively 
unobscured transient f or SN2008S Prieto et al.1 ( 20081) : 
Botticella et al. ( 2009!) illustrated this possibility, fol¬ 
lowed by the case of SN2012aw (Section !?.1.21) . The up¬ 
per limits suffer much more uncertainty than the direct 
detections as in the latter cases the extinction can be 
constrained by SED fits, at least for the cases with the 
largest wavelength coverage in the broad band colours. 
Nevertheless, under reasonable assumptions, the lim¬ 
its are a useful constraint. Approaches have been to 
adopt the best estimate of reddening toward the SN, to 
adopt an arbitrary extinction to represent the CS M, or 


to mo del it from theoretical assumptions. ISmart t et al 


( 20091) added an additional Ay ~ 0.3 for those progen¬ 
itors which had only Milky Way foreground estimates 
(the value of 0.3 coming from comparisons between the 
SN sa mple and the red supergiant population of the 
LMC). iFraser ( 2011 ) reanalysed this sample, using a 
better method of model star SEDs (MARCs models) 
and adopting extra extinction of Ay ~ 0.5 for all (again 
from comparing the M31, LMC and SMC populations 
of red supergiants). To be conservative the val ues for 
t he lu minosity limits in Tabic]?] are quoted from IFraser! 
( 201ll) with this higher, ad hoc (but physically plau¬ 
sible) extra extinction of Ay ~ 0.5. The results of the 
theoretical modelling of stel l ar mass- l oss ra tes and dust 
production of Walmswel l fe Eldridg e (2012) will be dis¬ 
cussed further in Section !?)?! 

2.2 Type lib supernovae 

As discussed in the previous two sections, the pre¬ 
discovery images of type II-P SNe are often not deep 
enough to detect the progenitor stars. This is not sur¬ 
prising since for distance moduli of /r = 30 — 32 (10- 
25 Mpc) and the typical depth of HST images (hub ~ 
25 m ), the lowest luminosity red supergiants would go 
undetected. If the broad picture of core-collapse, neu¬ 
tron star formation and successful SNe occurring in 
stars above an initial mass of ~8M 0 is true, then the 
initial mass function must dictate that the bulk of CC- 
SNe we detect are from the lower mass region. The 
lib SNe are those which begin by resembling type II 
SNe with unmistakable Hi Balmer l ines an d ev olve t o 
have H-weak and He-strong spectra ( Fihjajaenko 199ljb 


the prototype of which is SN1993J (iFilippenko et al 


1993 : Matheson et al. 2000). It’s remarkable that the 
three lib SNe in the 1999-2013 distance and time 
limited survey which have high quality imaging all 
have relatively bright progenitors (or progenitor sys¬ 
tems) detected at the SN po sition. The de t ection s 
of the progenitors of SN2008ax ICrockett et al. ( 20081) : 
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Table 2 Upper limits for the luminosity of type II and type IIP SNe. Most of the limits quoted here were re-calculated by Fraser 
(2011), based on the original published photometry limits and with an additional, but adhoc extinction of Ay ~ 0.5 mag applied as a 
conservative estimate. The values for the initial mass of star which ends it life with that luminosity is given in the last two columns. 
As in TablefTI these are for the STARS and Geneva models M(S,G) and the KEPLER models M( K), as discussed in Section [3TI 


SN 

log L 

Reference 

M(S,G) 

M(K) 

SN2006ov 

< 4.7 

From 

(1) ; <4.4 from (2) 

< 10 

< 12 

SN2004dg 

< 4.7 

From 

(1) ; <4.5 from (2) 

< 10 

< 12 

SN2001du 

< 4.7 

From 

(1) ; <4.7 from (2) 

< 10 

< 12 

SN2006bc 

<4.9 

From 

(1); <4.4 from (2) 

< 13 

< 14 

SN2007aa 

< 4.9 

From 

(1) ; <4.6 from (2) 

< 13 

< 14 

SN1999gi 

< 4.9 

From 

(1); <4.6 from (2) 

< 13 

< 14 

SN2009N 

< 4.9 

From 

(1) 

< 13 

< 14 

SN1999br 

< 5.0 

From 

(1) <4.9 from (2) 

< 15 

< 16 

SN1999em 

< 5.0 

From 

(2) 

< 15 

< 16 

SN2009ib 

< 5.0 

From 

(3); same result in (1) 

< 15 

< 16 

SN2002hh 

< 5.0 

From 

(1); <5.0 from (2) 

< 15 

< 16 

SN2009H 

< 5.1 

From 

(1) 

< 16 

< 18 

SN1999an 

< 5.2 

From 

(1); <5.2 from (2) 

< 17 

< 20 


m lFrased (1201 lh ; (2 lSmartt et al.1 (12009ft : (3 lTakats et all (120151) 


Table 3 Type lib detections 


SN 

T e s 

err 

logL 

err 

Comments 

SN2008ax 

3.95 

0.2 

5.1 

0.2 

(1) ; See discussion in Sectl2.2l 

SN2011dh 

3.78 

0.01 

4.9 

0.2 

From (2). Similar results in (3) 

SN2013df 

3.62 

0.01 

4.94 

0.1 

From (4) 


(T lCrockett et al. ( 2008! ); (2)Maund ct al. (20111); (3Wan PTet al. (201l|); (4fVan Dvk ct al. ( 2014 ): 


SN2011dh (Maundetal. 2011 


and SN2013df (Van Dvk et al 


Van Dyk et al. 1201 1); 


201 4) are summarised 


in Tabled The spectral energy distribution of the pro¬ 
genitor of SN2011dh is very well fit with a spectrum 
of a single star with T e s = 6000 K and logL = 4.9 Lq 
as shown i n Figured The prog enitor of SN2013df is 
cooler, with l Van Dvk et ~ah ( 20141) finding T e g = 4250 K 
and remarking that the detecte d prog eni tor appear s 


quite similar to that of SN1993J ( Alderin g et al 


ppear s 

Il994h . 


The progenitor of SN1993J was postulated, and shown 
to be, a binary system of a K-type supergiant and 
hotter comp anion. It too had a luminosity of log L = 


5.2 ± 0.3Lf7i (Van Dvk et all 12002: Maund et all 2004: 


Fox et ah 20141). 

The colours of the progenitor of SN2008ax could 
not be satisfac torily fit with a single stellar SED 
(Cr ocket t et al. 120081 ) or with a plausible binary which 


led the authors to suggest it might be a single WR-type 
star or a system which has contaminated flux from a 
physically associated cluster or nearby stars. Crockett et 
al. found the absolute magnitude of the star coincident 
with the SN position to be Mf 606W = —7.4 ± 0.3, and 
if we assume a bolometric correction of —0.6 (to cover 
possible supergiant stars from late B-type to early K- 
type) then this implies a stellar logL ~ 5.1Lq but with 
quite large uncertainties due to the difficulty in fitting 


a single SED. Groh et al. (2013b) suggest that the pro¬ 
genitor could have been a rotating star of around 2OM 0 
that ends its nuclear burning life as a “compact LBV” 
with a luminosity of around logL = 5.2 — 5.3 dex. Nev¬ 
ertheless, taken together the lib progenitors tend to 
be brighter, more often detected, and of systematically 
higher luminosity than the bulk of the type II-P SNe. 

The model of lib SNe coming from extended 
supergiants which have had mass stripped due 
to a binary companion was supported initially 
by the d etectio n _of a prog en itor comp anio n in 
SN1993J ( PodsiadlQwskiel_aL 1993; Nomoto_etjd 


199 3: Wooslev et al . 1994 iMaund et al . 2004; Ifox et al 


201 4) and recently also for SN2011dh At the time of dis¬ 
covery of SN2011dh, a yellow star dominated the SED 
of the progenitor system and the community had split 
views on whether the yellow supergia nt or i ts p u tative 
comp act companion had exploded (e.g. lSoderberg et al.1 
2 0121) . Th e disappea rance of the bright yellow star 
(IVan Dvk et al.ll2013l) and the de tection of a UV source 
remaining at the SN position bv iFolatelli et al.1 ( 2014 ) 
is strong evidence that it was indeed the yellow super¬ 
giant that exploded and tha t it had a more comp act 
companion (see the model in Benvenuto et al.ll2013 ). 

It remains to be seen if the existing, and follow-up 
data for SN2013df and SN2008ax also can be explained 
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5.5 5.0 4.5 4.0 3.5 

°9 T eff 


Figure 3. The positions of the detected progenitors and upper limits to the type II SNe as discussed in Section[2] The stellar evolutionary 
tracks are from lEldridge fc Tou t (2004). The possible positions of the progenito r star of PTF13bvn is ma rked with two symbol s, joined 
by the dotted line. These show the two positions of the progenitors proposed bv lBersten et all d2014f) an d El dridge et alj (120151) in their 
binary models. The position of the progenitor of SN2009ip is shown as the magenta symbol, as estimated from the faintest magnitude 
the LBV star was found at (see Section 2.4 for more details). The 14 Ibc progenitors with no detections are not quantitatively marked 
here. If they were WR stars then one would expect to find them around the blue shaded area (although the box position is illustrative 
as some models predict progenitors outside this locus, e.g. iGroh et al.ll2013a| [ch There are 30 progenitors below log L = 5.1, and only 
one (SN2009ip) above, if indeed SN2009ip is a genuine core-collapse supernova. 


with binary models. However it is certain that the pro¬ 
genitors all have luminosities around logL ~ 5.OL 0 , al¬ 
beit with the assumption that the bolometric correc¬ 
tion for the SN2008ax progenitor is not unusually high. 
The stellar evolutionary models employed in the stud¬ 
ies discussed in this section to explain the progenitor 
luminosities have ZAMS masses in the range 13-17M 0 . 
It would appear that lib progenitors are consistently 
among the highest luminosity (and hence highest mass) 
progenitors so far detected. There are few other upper 
limits, and with four cases of nearby lib SNe all with 
bright, detected progenitors it seems reasonable to con¬ 
clude that they are not produced by lower mass systems 
in the 8-12M 0 range. The physical reason for this is not 
determined and it will be interesting in the future to test 


quantitative models. One may speculate that it could be 
due to higher mass stars forming closer binaries more 
easily or brighter progenitors being more extended and 
hence their envelopes more efficiently undergoing Roche 
lobe overflow when they reach the LI point. 

A spectrum of the light echo of Cas A showed th is 


Milky Way supernova to be of type lib (|Krause et al 


2008 : Rest et al. 2008). There is no sign of an obvi¬ 


ous companion star in the centre of the Cas A rem¬ 
nant, leading to the speculation that this was a high 
mass single star that lost most of its hydrogen en¬ 
velope through winds. However the total amount of 
gas either ejected or residing in the CSM, and there¬ 
fore recently ejected by the progenitor is not partic¬ 
ularly high. The mass of the shocked and unshocked 
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ejecta is estimated at 3 — 4Mq ( Hwang fc Laminj2012 ) 
and the mass in the CSM which the blastwave has 
swept up is around 8 M^ ( Hwang fc Lamin3 20091 : 
Patnaude fc Fesenl 20091 ) . Adding in ~ 2M 0 for a neu¬ 


tron star remnant, results in a total progenitor mass 
estimate of about 15 Mq. The internal structure seen 
in the ejecta originate from turb ulent mixing processes 
and plumes of 56 Ni-rich material (Milisavl i evic fc Fesen 


20151) . These detailed constraints are in reasonable 


agreement with the the results from the three lib SNe 
with progenitors discussed above - originating from 
stars with initial masses in the range ~ 15 ± 3M 0 . How¬ 
ever the puzzle remains for Cas A - if the star had no 
companion, and the initial mass was only ~15M 0 , then 
how did it lose its envelope to become a lib SN? One 
might speculate it was a high mass, obscured star, al¬ 
though it seems the exti ncti o n towar d Cas A is no more 
than Ay — 6 — 7 mag ( Eriksen et al.l l2009h and this is 
integrated along the Galactic line of sight rather than 
being local to the SN and the progenitor. 


2.3 Type Ibc supernovae 


There are two possible models for the progenitors of 
type Ibc supernovae. The first are Wolf-Rayet stars 
which are evolved, single massive stars that have 
lost their hydrogen envelopes (in the case of WN 
stars) and also helium layers (in the case of WC and 
WO stars) primari l y through radiatively driven wi nds, 
( Kudritzki fe Puls 2000 : Vink fe de Koter l2005h or 
episodic mass-loss ( Smith fc Owockil 2006 ) or through 
rapid rotation and ch emically homogenous evolution 
(]Yoon fe Langenl2005l) . The observed Wolf-Rayet stars 
in our galaxy and the Magellanic clouds are observa- 
tionally constrained to come from stars (either single 
or binary) above about 25-3OM 0 from considerations 
of the turn off masses in stellar clusters (e.g. see re¬ 
view of lCrowther 200 7). From a theoretical standpoint, 
the stellar winds of single stars below this mass are 
not strong enough alone to cause sufficient mass loss 
to expose the He and CO cores. These stars are pre¬ 
dominantly found in young stellar clusters and OB as¬ 
sociations within the Local Group and turn-off masses 
for coeval stellar populations have provided estimates 


of mass and age (e.g. 1 Massey et a l. 


19951). The second 


alternative is lower mass stars in binary systems with 
i nitial masses lower than that required for single stars 
( Paczviis klll 1967: Vanbeve ren et al.lll998t Nomoto et, al 


19951: Podsiadlowski et, al. 1992 : lEldridge et,~ 2008 
Y oon et al. 20101) Theoretical stellar populations in¬ 
cluding binaries with reasonable distributions of stellar 
masses and separations can produce large number frac¬ 
tions of He stars with c ore masses large enou gh to un¬ 
dergo core-collapse (e.g. Eldridee et al. I20131) . However 
one problem with these models is that the systems they 



nitial Mass M/M 0 


Figure 4. A comparison of the final pre-SN luminosity as a func¬ 
tion of stellar mass, from three sets of stel lar evolution models, 
discussed in the text ll.Terkstrand et al.l2014l ). The STARS and ro¬ 
tating Geneva models have very similar core masses throughout 
the mass range, whereas the KEPLER models are 0.1 — 0.2 dex 
less luminous at the same stellar mass. This would corresponds to 
difference in estimated progenitor mass of 2-3 Mq if the final lu¬ 
minosity is used as an initial mass tracer. The progenitor limis for 
the ’’DUSTY” model are those for SN2012aw, from and are dis¬ 
cussed further in Section f3.5l Reproduced from Figure 5, in “The 
nebular spectra of SN 2012aw and constraints on stellar nucle¬ 
osynthesis from oxygen emission lines”, A. Jerkstrand et al, 2014, 
Monthly Notices of the Royal Astronomical Society, 439, 3649. 


predict are not observationally identified in the Milky 
Way. 

As discussed in Sectional t h e volu me and time lim¬ 
ited survey of Eldridge et al. ( 201 3) presented an ex¬ 
tensive search for the progenitors of all type Ibc SNe 
in pre-discovery images (from 1998 to 2012.25). This 
paper described the compilation of literature data, and 
new analysis of images of the progenitors of a total of 
12 type Ibc SNe. There are no detections of progeni¬ 
tors, or progenitor systems in any of these. The deepest 
limits in the typical BVR- bands were between —4 and 
—5 and, the authors compared these limits with the 
observed magnitudes of WR stars in the Large Magel¬ 
lanic Cloud. As WR stars, or stripped He stars in bi¬ 
nary systems, have diverse temperatures and radii (and 
hence diverse colours, SEDs and bolometric corrections) 
the quantitative fitting methods which have been suc¬ 
cessful for type II progenitors cannot be applied to the 
Ibc limits without major uncertainties in the results. 
Hence a more empirical approach was applied which 
simply postulated that the observed WR star popula¬ 
tion of the LMC are plausible progenitor systems of 
Ibc SNe and then determined what was the statistical 
uncertainty that no detections were made. Since then, 
there have been two further Ibc SNe within the distance 
limit, also with no detections of progenitors: SN2013dk 
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( Elias-Rosa et al. 2013) and SN2012cw ( Graur fe Mao d 
2012 ). Eldridge et al.1 ( 20131 ) originally estimated the 
probability of not finding a progenitor system detec¬ 
tion, if the LMC WR population are direct progenitors 
of Ibc SNe to be 16%, and with the inclusion of the 
limits of SN2013dk and SN2012cw this drops to ~12%. 
This is an interesting constraint but not in itself strong 
enough to rule out WR stars being the progenitors of 
the normal Ibc population we see in the local Universe. 
Several other authors have argued that the temperature 
and luminosity (which of course dictate optical and NIR 
fluxes) of the observed WR population are not applica¬ 
ble to WR stars at the point of collapse and that signif¬ 
icant evolution to higher temperatures and fai nter opti¬ 
cal magnitudes oc cur in their models ( Yoon et al. 20121: 
Groh et al. 2013d) . These of course are dependent on 
the mass loss rates and the time to core-collapse from 
where we observe them now. It is also possible that the 
extinction toward these progenitors, or CSM extinction, 
is systematically underestimated as discussed above for 
the type II-P progenitors and limits (in Section !?. 1.31) . 

There are two other strong constraints on the progen¬ 
itors of type Ibc SNe. The first is their relatively high 
rate which has been explained with binary population 
models (for example as disc ussed in lPodsiadlowsldjdjd 


1992;: Eldridge et al. 1200 8). More recentlv lSmith et*aT 


(|2011al) pointed out that the the SN rates from the 
volume limited Lick Observatory Supernova Search 
(LOSS) imply a rate (26% of all CC) that is a fac¬ 
tor of ~2 too high for massive, single WR stars to be 
the sole producers just from simple IMF arguments. An 
IMF of single stars would produce a WR fraction (and 
hence a Ibc fraction) of 14-18% of all massive stars be¬ 
tween 8-100 M^ if the minimum ma ss to produce a WR 
is 25-3OM 0 . lEldridge et al.1 ( 20131 ) built on this with 
the Binary Population and Spectral Synthesis (BPASS) 
code to calculate rates of progenitors in stellar popu¬ 
lation models including binaries. The rate of Ibc SNe 
within 28Mpc (during 1998-2012.25) is also 26%, and 
Eldridge et al. ( 20131 ) find that the rates can only be 
produced if a mixed population of single stars and bi¬ 
naries are used with one single star for every binary 
system and the mass ratio and initial separations are 
set such that approximately two thirds of the binaries 
interact. This valu e is co mpatible with the observed re¬ 
sult of Sana et al. ( 201 2) who showed that this fraction 
of massive binaries are likely to interact. These multi¬ 
ple and consistent observational constraints from rela¬ 
tive SN rates and the binary fraction of massive stars 
are rather compelling arguments that binary systems 
are likely to produce the bulk, if not all of the Ibc SN 
progenitors. lEldridge et al.1 1 2013 ) points out that the 
reason that Galactic analogues may be difficult to iden¬ 
tify is that stars in binaries with initial masses below 
20Mq tend to retain a low-mass hydrogen envelope un¬ 


til quite close to core-collapse, ~ 10 4 yrs. Hence they are 
not easily detectable as H-free compact helium stars. 

The second additional argument for Ibc SNe coming 
from interacting binaries is that the ejecta masses cal¬ 
culated from ligh tcurve modelling are typically in the 
range 1-4M 0 (see I Valenti et al. 2008; Drout etal]l2011 . 


the results of which are summarised in Eldridge et al. 
2013). These final masses are somewhat low compared 
to the estimated masses of WR stars in the Local Group, 
which are typically in the range 8-2OM 0 (from binary 
orbits, Fig. 4 of lGrowtherl 2007 ). The binary models of 
Eldridge et al. (2013) produce helium stars that would 
have typical ejecta masses of 4.2 ± 2.4M^, in agre ement 
w ith those measur ed (also see Lvman et al. 20141). 


Anderson et al. ( 2012 ) and Anderson fe Jamesl 


(120081) have matched the positions of nearby SNe to 
the Ha emission produced in starforming regions. The 
SNe Ibc in their sample are more closely associated 
with the Ha emission than type II SNe in general. 
They interpret this as the SNe Ibc (and SN Ic in 
particular) arising from younger re gions and hence 
coming from higher mass progenitors. iGrowther ( 201 3) 
compared a much closer sample of CCSNe (<15Mpc) 
with higher resolution images and also concluded that 
the association of Ibc with Ha regions was stronger 
than for type II SNe. However the diagnostic power 
of this association in quantitatively determining the 
mass ranges was shown to be difficult to interpret. 
The Ibc progenitors do appear to come from younger 
populations than the bulk of the type II progenitors. 
However the latter are dominated by stars in the 
8-12M 0 regime, and a quantitative age and mass 
estimate for the Ibc progenitors has not yet been 
possible. Although they may on the whole be from 
more massive progenitors, the data can’t at present 
distinguish between WR stars and lower mass (say 
12-2OM 0 ) binaries. 


2.3.1 PTF13bvn 

Just after the publication of the Eldridge et al. (2013) 
summary, the Palomar Transient Factory discovered 
a nearby lb SN in NGC5806 (d = 22Mpc) an d set a 


very restrictive time on the explosion epoch (Cao et al. 


20131) . The non-detections indicate that the first detec¬ 
tion was within 24hrs of the shock breakout. The HST 
pre-explosion images of the site indi cate a blue star at 
the position of the SN. iGao et aD ( 2013i ) noted that 
their alignment of a high resolution ground-based image 
of the SN was close to, but not formally within the lcr 
error circle of the alignment calculation. However one 
cannot confidently reject a physical association unless 
there is a 3a difference and the a is co rrectly defined 


and m easured. Subsequent alignment in Eldr idge et al 


(12015 ) with HST imaging again found a simil ar re - 
sult but supported the suggestion bv ICao et al. ( 2013 ) 
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that this was likely the progenitor system and the 
first detection of a stellar progenitor of a Ibc SN. 
Cao et al. ( 2013 1 and lGroh et al. ( 2013all proposed that 
the broad band magnitudes of the source in the HST 
images were similar to massive WR stars and a stel¬ 
lar model of initially M z ams — 3OM 0 evolved into a 
WN star with model parameters logL/L odot = 5.6 and 
T e g = 46kK (producing a final CO star mass of 1OM 0 ). 
A reassess ment of the HST ma gnitudes of the pro¬ 
genitor led lEldridee et al.1 ( 20151 ) to argue that a bi¬ 
nary system could explain the updated stellar mag¬ 
nitudes better. They found brighter magnitudes by 
about 0.5 m ; -Mf 435 W = — 6.0 ± 0.4, -MF 555 W = —6.0 ± 
0.4, Mf 814 w = —5.9 ± 0.4. These updated magni t udes 


do not rule out the single WR models of Groh et al 


(] 2013a) . but the model binary systems reproduce the 
measur ements q u ite w ell. 

Both Bersten et al. ( 2014 ) and Eldridee et al.1 ( 2015 ) 
presented binary models as alternatives which have ini¬ 
tial masses of either 20 + 19M 0 or 10 + 8M 0 compo¬ 
nents. These produce a stripped He-core with a final 
mass of Mfi na i = 3 — 4M 0 and the combined flux from 
the binary systems can reproduce the progenitor magni¬ 
tudes. The ejecta mass would th en be of order 2M^ in 
agreement with the estimate of iFremling et al.l (J2014I ) 
from lightcurve modelling. While a single, massive WR 
star cannot be definitively ruled out, the progenitor 
magnitudes and ejecta mass estimates are quite well 
reproduced with a binary system (see Section liPl for a 

position on the Hertzsprung Russell Di agram). _ 

In a recent development in this area, I Gal-Y am et al 


( 20141 ) took a spectrum of SN2013cu within the first 
day after explosion showing that it has spectral fea¬ 
tures similar to those seen in WR stars (specifically 
Hi, He 1 , Hen, Civ and Niv seen in WN type stars). 
While this could be inter preted as evidence for a di¬ 
rect WR progenitor star, Gal-Yam et al. ( 20141) note 
that the inferred wind density and mass-loss rate are 
rath er extre me when compared to typical WR stars. 
Groh ( 20141) applied a more detailed analysis code to 
suggest that a WR progenitor star, with a classical ra- 
diatively driven wind, is not likely and the precursor was 
possibly an LBV-type or yellow hypergiant with a re¬ 
cent eruptive phase. SN2013cu was a type lib and the 
mass-loss event that produced the circumstellar mass 
illuminated in the wind could have plausibly occurred 
during mass-transfer in a binary or in an eruptive phase. 
The remarkable spectrum is similar to Galactic WR star 
spectra, although the emitting surface is likely a factor 
~10 larger than the typical WR stellar radii and the 
the narrow cores of the emission lines seen in SN2013cu 


are no t ty pical o f fast W R wind lines. As ISmith et al 


(l2015l) and[G roll! d2014l) point out, the spectrum of the 
early SN is not a direct indicator of the spectral type of 
the progenitor star before explosion, but it is a powerful 
measure of the recent mass-loss history of the progen¬ 


itor. It shows that the CSM material within ~ 10 AU 
is similar in H and He composition to WN star winds, 
rather t han the sta r definitely being a WR-type pro- 


2015). 


genit or (|Grohl 12014 iShivvers et al.l 12014 ISmith et al 


2.4 Type Iln progenitors : LBVs and 
SN2009ip 

There are three SNe which have been proposed as core¬ 
collapse SNe of type Iln with either identified high lumi¬ 
nosity stellar progenitors detected or reasonably strong 
arguments in favour of high-mass progenitors (other 
work has suggested LBV-like progenitors from modula- 
tions in radio lightcurves or in optical sp ectral features 
Kotak fc Vinkl[200fl iTrundle et al.ll2008l ). The first dis¬ 


covered was the progenitor of SN2005gl which was a 
bright My ~ —10 progenitor de tected at the surpris¬ 
ingly large distance of 60 Mpc (IGal-Yam et al.1 12007 : 
iGal-Yam fe Leonard 20091) 

The type Iln SN2010jl has HST pre-explosion imag¬ 


i ng of its host galaxy (at 50Mpc) and Smith et al 


(2011b) identified the site of the explosion. At this dis¬ 
tance the relatively crowded region of the explosion site 
prevents a confident detection of a resolved point source, 
but the blue flux coincident with the SN is likely to be 
either a ver y young cluster or a single massive star. 


Smith et, all (1201 lb ) argued that either of these two 


cases suggest that the progenitor had an initial mass 
higher than Mzams > 3OM 0 . Both of these SNe oc- 
curred within the t ime l imit of the survey de fined by 
ISmartt et al. ( 20091) and Eldridee et al. ( 2013 ) but are 
much more distant than the 28 Mpc limit. 

The now infamous SN2009ip is a type Iln SN 
( Mauerhan et al. 20131 ) with a progenitor detected in 
outburst over many years. After the initial discov¬ 
ery of a non-t erminal bu t luminous (M ~ —14) out¬ 
burst in 2009 ( Maza et al. 20091) at which point it was 
named SN2009ip, a luminous progenitor star {My = 
— 10) wa s detected in HST archival i magin g taken 
in 1999 ( Smith et al. 2010t Folev et al. 201 lh . Subse- 
quent monitor ing over the next three years, mostly by 


Pastorello et al. (120131 ) illustrated that the explosion in 


SN2009ip was part of a regular, but temporally spo¬ 
radic, series of outbursts of a massive star which photo¬ 
metrically shared similarities with the Luminous Blue 
Variables (LBVs) of the Local Group. Indeed it is not 
even clear if the detection in the HST imaging is the 
star in outburst or quiescence and hence attaching a 
bolometric luminosity to any of the data points is not 
easy to interpret with hydrostatic stellar evolutionary 
models. It may be that the luminosity is powered at all 
phases, including that in the archival HST images, by 
interaction (and the same could be true for SN2005gl). 

In 2012, the star had a 30 day outburst that 
peaked at a luminosity similar to the brightest known 
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LBV giant outbursts (the“2012a” event ; My ~ —15, 
logL = 41.5ergs -1 ), which was followed by a 10 day 
rise to SN luminosities (the“201 2b” event; My ^ —18 . 
log A = 43ergs - G (monitoredby Pastorello et al.llioi.St 
Mauerhan et al. 20131 : IPrieto et al. 20131). The cause 


and nature of the two are still debated. Smith et al 


( 20141) fervently argue that core-collapse must have oc¬ 
curred to cause the bright 2012b peak, and indeed many 
observational characteristics of this explosion are simi¬ 
lar or ident ical to type I l n sup ernovae. 


Mauerhan et al 


(2013) and S mith et al. (2014) propose that the fainter 
“2012a” event was a non-terminal giant eruption, that 
core-collapse caused gas to be accelerated to veloci¬ 
ties of > 10,000 km s -1 and the kinetic energy of this 
gas is efficiently converted to radiative energy in col¬ 
liding with a dense, but p ossibly asymmet r ic CS M to 
produce a bright SN . Both IPastorello et al.l (|2013l) and 
iFraser et al. ( 2013al) argued that it is still possible that 
the last SN-like event in the history of SN2009ip was 
not due to core-collapse, but caused by the collisions 
of massive shells ejected from the star with dense a 
CSM. T he energetics are plausible, although are dis- 


Graham et al- ( 2014 ): 


puted bylSmith et al. (2014 ). Other papers o n the topic 


Margutti et al.l ( 2014 ) have also 


presented extensive data and analysis but none as yet 


have definitively proven the true nature. Ofek et al 


( 20141) have suggested that at least 50% of type Iln 
SNe have outbursts within the last few years of their ex¬ 
plosion. This seems plausible given the Palonrar Tran¬ 
sient Factory (PTF) discovery r ate, the other discov¬ 
eries of pre-explosio n outbursts (IPastorello et, al.ll2007 : 
Fraser et al.l l2013bl) and the theoreti cal ideas of hy¬ 


drody namic instabilities discussed in ISmith fc Arnett 
(20141). Whatever mechanism powers the bulk of Iln 
SNe it appears clear that eruptions in the few years be¬ 
fore explosion are much more frequent that previously 
thought. 

For the purposes of this review, the interesting ques¬ 
tion concerns the possibility that very massive stars 
produce successful SNe through the core-collapse mech¬ 
anism. The faintest magnitude that the progenitor star 
of SN2009ip was observed at was My = —10, and with a 
bolonretric correction appropriate for a blue supergiant 
(—1.0 <; BC <j —0.5), this equates to logL/L 0 = 6.0 - 
6.3 dex, or an initial mass of 6O-8OM 0 . The discussion 
that follows, will consider the possibility that SN2009ip 
has undergone a core-collapse and successful explosion. 

Although it has been reasonable to link supernovae 
of types Iln to LBVs or very massive stars that have 
had violent recent mass-loss history, an alt ernative ex¬ 
planat ion was recently put forward by Mackev et al.1 
(12014 ). In this model, a red supergiant with a fairly 
standard mass-loss rate (such as Betelgeuse with M = 
1.2 x lO -6 M 0 yr -1 ) produces a static shell of circum- 
stellar material through pressure from photoionization 


by radiation from external sources. The steady wind 
gets confined into a static photoionization-confined shell 
which could contain up to 35 per cent of all mass lost 
during the red supergiant phase. This gas shell is con¬ 
fined close to the star, giving a possibility of a type Iln 
supernova resulting from a moderate mass red super¬ 
giant with standard mass-loss. 


3 DISCUSSION 

3.1 The luminosity distribution of 
progenitors 


A summary of the detected progenitors and their lim¬ 
its are illustrated in a Hertzsprung Russell Diagram in 
Figure[3l along wit h a set of theore t ical s tellar evolu¬ 
tionary models of lEldridge fc Tout! ( 2004 ). This rep¬ 
resents a complete summary of the progenitors de- 
te cted (and limits set) withi n the time and volume limit 
of Smartt et al. ( 20091) and Eldridee et al.l ( 20131) (ex¬ 
tended to all data published up to end of 2013) and 
hence the relative numbers and distribution of the stars 
in this plot are meaningful. The striking result is that 
there is a deficit of progenitors above an estimated lumi¬ 
nosity of log A/L 0 ~ 5.1. There are 30 progenitors (type 
II detections, II limits, lib and lb detections) below 
this luminosity limit. The final luminosity of a super¬ 
nova progenitor is determined by the He core mass and 
luminosity, and this luminosity equates to a slightly dif¬ 
ferent mass depending on the stellar evolutionary mod¬ 
els used. This log I//Lm~ 5.1 equates to 16Mnfor the 
STARS models of lEldridge fe Tout (12004) : 16M^ fo r 
the rotating Geneva models of Hirsc hi et al.1 (120041) ; 
18M 0 for the IWooslev fc Heger (120071 ) (see Figured]) 
and other stellar evolutionary codes produce model pro¬ 
genitors that are broadly in this range. 

For a Salpeter IMF with a = —2.35, 70% of stars be¬ 
tween 8-lOOM 0 are in the 8-18M 0 range. Hence with 
this number of stars (30) lying below the mass limit 
of 18M 0 one would expect to have found 13 higher 
mass progenitors with log A/L 0 > 5.1 dex and Mzams > 
18M e 3 There is possibly one (SN2009ip; Section ETfl) . 
but the Poissonian probability of finding between 0 and 
1, if the expectation value is 13 would be p = 3 x 10 -5 . 

The number 13 is a very close match to the number 
of Ibc progenitors for which we have upper limits and 
one detection. One could argue that the bulk of stars 
above 18M 0 evolve into WR stars and evade detection 
due to them being too hot and faint at the point of 
core-collapse. And the numbers would match (almost 
too well). This is not strictly ruled out by any of the 
observational constraints to date but it appears to be 


1 Or alternatively one could propose that these 30 objects are the 
total sample and actually do contain the full range of masses. 
This would mean around 9 stars of this sample would need to 
have had their luminosities significantly underestimated. 
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unlikely. The three constraints of having no detected 
progenitors, the high relative rate of Ibc SNe and the 
low ejecta masses estimated from liglitcurve fitting (pre¬ 
sented in Section l2T3l) all suggest that the binary channel 
dominates the production of type Ibc SNe. 

The missing high mass stars clearly exist in sub¬ 
stantial numbers in starforming regions and we observe 
them following a Salpeter (or similar) initial mass func¬ 
tion. Therefore the question is, why are they not de¬ 
tected as SN progenitors ? It does not appear likely 
that they are all WR stars that have evaded detec¬ 
tion, given the arguments above. One argument put 
forward is a bias in detecting high mass stellar pro¬ 
genitors as they may exist in denser, younger star¬ 
forming regions - dense clusters or associations. Thus 
the surface brightness of these regions could make re¬ 
solving individual stars more difficult. However there 
are only three known SNe within the distance and 
time limit that fall on compact clusters. Two of them 
SN2004am and SN2004dj are relatively old clusters 
and very unlikely t o host stars of Mzamr > 2 0M^, 
( Mattila et al. 120131 : iMafz- Anellaniz et al. 120041) and 
the other is SN1999ev which is also unlikely to be a 
very young cluster ( Maund et al.l 2014bl) . It is also un¬ 
likely that these missing high mass progenitors stars 
produce normal luminosity SNe that are systematically 
missed by nearby surveys because of surface bright¬ 
ness effects - since SN2004am was a relatively faint 
and reddened SN, rec overed on top of a bright super 
star cluster in M82 ( Singer et al. 20041 iMattila et al 


20131) . And S N2004dj was also easily detected in ama¬ 
teur surveys (jNakano fc Itagakill2004l) . It would appear 
that the SN population we observe in the local Universe 
are not being produced, in large proportions, by stars 
with logL/L 0 > 5.1 (M Z ams ^ 18M 0 ). 


3.2 The mass function of progenitors 

Smartt et al. ( 20091) calculated the mass function of 
progenitors of type II SNe, assuming that there are 
no major biases in selecting the SNe with progen¬ 
itor information from the broader local SN popula¬ 
tion. This paper used a consistent approach, employ¬ 
ing the STARS models to attach an initial stellar mass 
to the luminosities and luminosity limits determined. 
They defined the “red supergiant problem” which ar¬ 
ticulates the fact that while red supergiants are com¬ 
monly found as direct progenitors of II-P SNe, there 
were no progenitors found above an estimated initial 
m ass of around 17M 0 , o r logi/L 0 ~ 5.1. As discussed 
Levesque et al.| ( 20091) . red supergiants are found up 


to luminosities of logi/L 0 ~ 5.5. 

It is timely to review this in light of the discoveries 
made since 2009 and the increase in statistics now avail¬ 
able and updated values using different methodologies. 
The luminosities compiled in this paper for detections 


MF=- 1.35, Salpeter 



7 8 9 10 11 

Minimum Initial Mass / M 0 


Figure 5. The maximum likelihood of the minimum and max¬ 
imum initial masses of the type II progenitor distribution, as¬ 
suming the stars follow a Salpeter IMF mass function. Originally 
calculated in in Smartt et al. (2009), and reproduced here with 
the updated and extended masses in this review. The dashed lines 
show the confidence contours (68, 90 and 95%) for the detections 
only and the solid lines show the confidence contours for the 
detections and upper limits combined. The star symbol marks 
the best fit, as described in Section[3]2] of m m ; n = 9.5y)('~' and 
mmax = 16.5^2 5 - This is for the masses from the STARS and 
Geneva rotating models, the values for the KEPLER masses are 
given in the text. 


and non-detections are used with three sets of mod¬ 
els to estimate initial masses of the progenitor stars. 
The end points of three sets of models are plotted in 
Figure^] showing the initial mass - final luminosity re - 
lation for the STARS models [Eldridge fe Tout (I2004I) . 
the rotating Geneva mo dels oflHirschi et al. (2004) and 
the KEPLER models of Wooslev fc Heger ( 2007 ). The 
STARS and rotating Geneva models have quite simi¬ 
lar final luminosities, hence we consider the results of 
these to be indistinguishable (at least at the level of ac¬ 
curacy possible for any individual SN progenitor). The 
KEPLER models are less luminous at their endpoints 
by between 0.1 — 0.2 dex which typically translat es into 
a difference in mass of 1-2M 0 . In ISmartt et al.l ( 20091 ) 
and other work that employed one set of stellar mod¬ 
els, the uncertainty in the theoretical models was taken 
into account by using luminosities at the end of the He 
burning. This was somewhat artificial and this review 
uses an alternative approach of taking the end points of 
three sets of evolutionary model grids. 
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Figure 6. The progenitor detections are marked with error bars 
(data from TablefTI and the limits are marked with arrows (data 
from Table[2). The lines are cumulative IMFs with different min¬ 
imum and maximum masses. 


The estimated initial masses of the progenitors are 
listed in TablesQ] and [2] The masses for the three low¬ 
est luminosity progenitors (SN2003gd, SN2005cs and 
SN2009md) are uncertain in the KEPLER models, 
as this code does not yet have published endpoints 
for stars in the 8 -IOMq regime. This is a problem¬ 
atic regime in that 2nd dredge up may occur, pushing 
the progenitors to higher luminosities th an t heir more 

Eldridge et all 


massive counterparts as discussed in 


( 2007 ): Smartt et al.1 ( 20021) and references there in. The 
quoted, but uncertain estimates from employing the 
KEPLER code are simply the extrapolated differences 
from the STARS estimates. 

Smartt et al.l ( 20091) , one can assume that the 


As 


progenitors come from a mass function of some slope 
and take a = —2.35 as a reliable standard. A maximum 
likelihood calculation can produce an estimate of the 
most likely lower mass and upper mass of the distribu¬ 
tion, assuming the masses follow a Salpeter function. 
For the masses estimated with the STARS (and rotat¬ 


ing Geneva) models, the values determined (using the 
same IDL routine as employed in Smartt et al. 2009) are 
a minimum mass for the distribution of m m ; n 


= 9.5l°' 5 


and a maximum mass of m ma x = 16-5_2' 5 where the 
errors are the 95% confidence limits (see Figure[5]). If 
we employ the KEPLER models, then the values are 
TOmin = an d a maximum mass of m max = 18.5+4 

(again with 95% confidence limits). These results are il¬ 
lustrated further in Figure[6] where the masses are plot¬ 
ted with a Salpeter IMF (cumulative frequency func¬ 
tion). The plots show that the mass distributions are 
comfortably reproduced with a standard IMF between 
the lower and upper mass limits from the maximum 
likelihood calculations, but they need to be truncated 
at the higher mass. If one allows the mass function to 
vary up to say 30Mq, then the mass distribution cannot 
be reproduced. This is the same basic result as shown in 
FigEl" the population of progenitors is missing the high 
mass end of the distribution, but this time the IMF is 
quantitatively considered. The maximum likelihood cal¬ 
culation is visualised in this cumulative frequency plot 
- given an IMF slope, the line fit should go through the 
error bars of the detections and not conflict with any of 
the upper limits. 

The lower mass lim it to prod u ce a c ore-collapse SN 


was estimated in Smartt et al. ( 20091 ) to be 


to n 


8.5+ 10 


_4 - from t he sa me method and the sample to that 
point. Smartt ( 20091) reviewed the limits from the max¬ 
imum masses of white dwarf progenitors, suggesting a 
convergence at TO m i n = 8 ± 1. The two values estimated 
here slightly higher: the value from the STARS models 
is TO m i n = 9.5l°' 5 (integer mass models evolved through 
C-burning down to that mass have been calculated) 
which is not significantly different to that estimated 
previously given the errors. The value from the KE¬ 
PLER models is higher again, at TO m in = lOl^'s- How¬ 
ever low mass models (7-lOM 0 ) are not available from 
KEPLER and the values in this luminosity range were 
estimated assuming the same differential in luminosity 
between KEPLER and STARS models exists between 
7-lOM 0 as at 11M 0 (FigurelU). This is uncertain and the 
lower mass from KEPLER should not be treated as a 
quantitative estimate : TO m i n is critically dependent on 
the mass estimates for the three lowest luminosity pro¬ 
genitors and if these are adjusted down by ~lM 0 then 
the value of m m i n = 9.5 would be reproduced. Some fur¬ 
ther quantitative modelling of stars in this interesting 
mass range is required to reproduce the stellar luminosi¬ 
ties and produce either a Fe-core collapse or O-Mg-Ne 
core that collapses through electron capture. Despite 
this uncertainty at the lower end, the existence of a 
high mass upper limit for type II SNe appears be secure. 
The value is model dependent of course, but the basic 
result is that type II progenitors are statistically lacking- 
above a logL/L 0 ~ 5.1 dex. The final model huninosi- 
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ties at this value are stars with 77i max = 16.5 (19M 0 
at 95% confidence) for the STARS and Geneva models 
and m ma x = IS.5 (21.5M 0 at 95% confidence) for the 
KEPLER models. 

While Smartt et al.l ( 2Q09h first discussed this as 
the “red supergiant problem”, the lack of detected high 
mass progenitors is now a broader issue for all SN types. 
This broader issue of a lack of high mass proge nitors 
generally was discussed in IKochanek et ah ( 2008h . who 
took all historical and literature limits to that date. The 
fact that there are now three type lib SNe in our sam¬ 
ple with progenitor detections (Table[3|) and these also 
have luminosities less than ~5.1 dex illustrates that the 
missing high mass star problem is relevant for all type 
II SNe (type II-P, II-L, lib). We address the issue of 
type Iln progenitors below. 

Alternative methods of probing the mass function of 
SN progenitor stars have been recently advan ced by 
Williams et al. ( 2014 ) and lJennings et ah ( 2014 ). These 
use similar methods of quantifying the stellar popula¬ 
tion around historical SNe in galaxies closer than 8 Mpc 
and around SN remnants in M31 and M33. They use 
high quality Hubble Space Telescope imaging and care¬ 
ful stellar photometric measurements to determine the 
luminosity and masses of stars within t he i m mediate 
vicinity of SNe (typically within 50pc). Jennings' et al 


( 20141 ) find that there is a lack of high mass stars close 
to the SN remnants in M31 and M33, suggesting ei¬ 
ther the highest mass stars do not produce SNe, or 
that SNR surveys are biased against finding objects in 
t he ve ry younge st (< lOMyr old) starforming regions. 
Williams et a l. ( 20141) also suggest that their results 
are compatible with progenitors all coming from masses 
M < 20Mq although the uncertainties do not rule out 
the possibility of no upper-mass cutoff. 

3.3 Possible explanations 

The reasons for these missing high mass progenitors are 
discussed as follows 

3.3.1 Dust formation and circumstellar extinction 
As discussed in Section ^. 1.31 the extinction toward the 
progenitors is often estimated from the extinction to¬ 
ward the SN itself, or the nearby stellar population. The 
former estimates may not be directly applicable since 
the circumstellar dust around the progenitor stars can 
be destroyed in explosions - as in the case of SN2012aw 
and SN2008 S. 

Walmswell fe Eldridgel ( 2012 ) calculated the dust 


that could be produced in red supergiant winds and 
the extra extinction that this w ould produce. The idea 
is well motivated and valid, but IKochanek et al.l (2012) 
showed that treating CSM extinction with a slab of 
IS M material is not phy sically consistent. As shown 
in Kochanek et al. ( 20121 ). the progenitor of SN2012aw 


was thought to be quite a high mass star but correct 
treatment of radiative transfer in a spherical dust shell 
reduces the progenitor luminosity limit while comfort¬ 
ably fitting the optical, NIR and MIR detections and 
limits. The major concern for this sample is that the 
objects with limits only are effectively unconstrained 
- one could propose that these are all high mass ob¬ 
jects enshrouded in dust and having no detection at any 
waveband does not allow meaningful constraints. How¬ 
ever there is another important constraint which comes 
fr om x-ray and ra dio ob servations of type II SNe. 


Chevalier et al. (120061 ) considered a sample of six type 


II-P SNe with x-ray and radio observations. They esti¬ 
mated the mass-loss densities in the progenitor star’s 
stellar wind from the thermal x-ray and radio syn¬ 
chrotron flux which originates when the fast moving SN 
ejecta interact with a pre-existing, lower velocity stel¬ 
lar wind from the progenitor star. The results for these 
six SNe (which are in the compilation presented here, 
either in TablcQ] or [2]) , were consistent with the mass- 
loss expected for pro genito r red s upergi ants in the mass 
range ~ 8 — 2OM 0 . IDwarkadas ( 20Til) has now com¬ 
piled all the available literatur e x-ray detect ions fo r type 
II-P SNe and used the Chevalier fc Franssonl ( 2003 ) 
treatment of the free-free x-ray luminosity from the 
shocks produced by the SN ejecta and the stellar winds. 
Red supergiant stars have observed mass-loss rates 
that span a wide range of values from < 10~ 6 M ( n yr~ 1 
up to ~ lO~ 4 M 0 yr _1 ( Mauron fc Jossellnll201ll ). The 
observed x-ray luminosity is L x _ ray oc ( M/v ) 2 (where 
M is the stellar mass loss rate and v is the stel¬ 
lar wind velocity; Chevalier & Fransson 2003). With 
a wind velocity of around ~10kms -1 , the x-ray lu¬ 
minosities imply mass-loss rates of less than M < 
lO -5 M 0 yr -1 . The mass-loss rates of red sup ergiants 
are co rrelated with stellar lumin osities and D warka das 
( 20141 ) then used the iMauron fc Josselin ( 2011 ) re¬ 
lations (M = 4.7 x 10 _6 (L/10 5 ) 1 ') to illustrate that 
type II-P SNe do not have high enough x-ray lumi¬ 
nosities to have had progenitors significantly greater 
than log£/ L(T)~ 5.2dex or 19M<^. Similar arguments 
were used in IKochanek et ah ( 20121 ) specifically for the 
case of SN2012aw to show that the steady wind of 
M <) 10 -5 5 to 1O~ 5 M 0 yr _1 was consistent with the 
x-ray and radio fluxes measured. The existence of ob¬ 
scuring CSM dust cannot exist without CSM gas, in the 
form of stellar winds, with typical dust masses around 
~ 0.5% that of the gas masse s (from the typical opac¬ 
ities in Kochanek et al.l 120121) . Therefore an obscured 
progenitor should be brighter in x-ray and radio t han 
the typical luminosities observed to date. IDwarkadas! 
( 20141 ) shows that the lib and IIL SNe have typically 


higher x-ray luminosities than t he I I-P, with val ues 
around ~ 10 39 5 ergs _1 (also see Dwarkadas fe Gruszkoi 
120121) . This is consistent with their higher luminosity 
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progenitors, in the range of log-L/L^ ~ 5.1 — 5.2dex. 
As remarked upon by Dwarkadas (2 0141 ) the general 
agreement between the two independent methods is re¬ 
markably consistent. While dust obscuration should still 
be considered carefully for any progenitor and limit, it 
appears the sample as a whole is not significantly bi¬ 
ased and there are no major discrepancies between the 
optical-NIR luminosity and mass estimates of progeni¬ 
tors and the SN x-ray and radio fluxes. 

It is interesting to note that it is not just steady wind 
mass-loss that has been proposed to produce dense CSM 
shells, an d pote ntially ext ra ext inction, around progen¬ 
itors. Shiode fc Quataert (2014) have argued that after 
core neon burning, internal gravity waves can transport 
a super-Eddington energy flux out into the stellar enve¬ 
lope and cause a mass ejection of around IMq of mate¬ 
rial. They suggest that this could occur preferentially in 
progenitors with Mzams ~2OM 0 within a few months 
to a decade of core-collapse, with the most intense mass- 
loss occurring closer to core-collapse. This timescale is 
quite similar to the typical time differences between SN 
explosions and dates of progenitor detections. The dates 
of the observations of progenitors, or limits (for the SNe 
in Tables 1, 2 and 3) have a range of 2 months to 17 
years before the SN discovery dates (which one can as¬ 
sume is close to the explosion date to within ~weeks at 
worst). The median and standard deviation of the time 
of data being taken before SN discovery are 65 and 54 
months. It is possible that these ejection episodes could 
create dust shells, however the gas masses would be in¬ 
consistent with the x-ray and radio fluxes, as discussed 
for the steady mass-loss scenarios. There is no obvious 
correlation between probability of discovery and date 
of data taken, and some of the discoveries have images 
taken within 5-9 months before explosion (e.g. 2003gd, 
2005cs, 2008bk). 

3.3.2 Biased luminosity estimates and evolution to 
the blue 

An argument could be made that the observed limit 
that we see for red supergiant progenitors of logT/L 0 ~ 
5.1 dex is not significantly lower, than the high¬ 
est luminosity red su pergiants of log L/L 0 ~ 5.45 dex 
( Levesane et al.ll2f)09t) and that uncertainties in the ex¬ 
tinction and luminosity estimates from small numbers 
of broad band fluxes are high enough to argue that this 
is not a discrepancy at all. Rat her, t his indicates broad 
and basic agreement (see Davies et al.1 1201.4 for a dis¬ 
cussion on careful model atmosphere techniques). How¬ 
ever the increase in statistics now available, and the 
statistical and quantitative deficit of high luminosity 
progenitors as calculated in Section EOl would suggest 
that this is not a satisfactory answer. 

More broadly, the question then is - if there is no 
real “red supergiant problem”, then where Hertzsprung 
Russell Diagram (HRD) do these Mzams )> 18M 0 stars 


evolve to before the point of core-collapse and what SN 
do they produce ? They must end up with an iron core 
and, within the bounds of current theory, they must 
collapse. It does not appear that they can all produce 
WR stars and Ibc SNe for the reasons discussed in Sec- 
tion l2.3l The search for progenitors is not biased against 
any particular SN type, and the one II-L progenitor 
(SN2009kr), and the three lib progenitors do not have 
luminosities higher than the observed limit of 5.1 dex. 
If high mass stars (M > 18M 0 ) evolve to produce visi¬ 
ble supernovae explosions, there is no obvious reason we 
should not detect both the supernovae and their progen¬ 
itors, unless the SNe themselves are intrinsically faint 
and have evaded either detection or recognition. 

The type Iln SNe have been suggested as a way out. 
High mass stars, with high mass-loss rates produce high 
enough CSM gas masses that the SNe are observed to be 
Iln types. This is reasonable and indeed corroborated 
by tw o detections of likely high mass progenitors for Il n 
SNe ( Gal-Yam fc Leonard! 120091 : ISmith et~ah 2011bll . 
but we have not systematically found eno ugh of them 


to acco unt for the missing mass range. Smith et al 


(j2011al) argue for quite a high relative rate of Iln SNe 
in the LOSS survey and that these would naturally 
have high mass progenitor stars. It seems almost cer¬ 
tain that SNe Iln do come from high mass progeni¬ 
tors (see Section HOt . Some have detected progenitors 
(as in SN2009ip and SN2005gl) and they tend to have 
the largest x-ray and radio fluxes that require mass-loss 
rates which can not b e realis tically produced by low to 
moderate mass stars ( Dwarkadasi 2014 , and references 
therein). Ho wever t he high relative rate estimated in 


Smith et al.1 (l2011aT) does not appe ar to hold at lower 
distances (see lEldridee et akll2013l for a discussion on 
this topic). There is no reason why we should not de¬ 
tect them and no reason why they should not appear 
in Figure[3] other than there are not enough high mass 
stars producing visible Iln SNe. 

Other work has also hinted that there could be some 
systematic bias underlying the mass estimates from the 
direct detection of progenitors. Hydrodynamic mod¬ 
elling of type II-P SN lightcurves and velocity measure¬ 
ments has suggested a significant discrepancy between 
the SN ejected masses and the stella r masses of the pro¬ 
genito rs. T he detailed modelli n g of lUtrobin fe Chugai 
( 20081) and lUtrobin fe Chuea'i ( 2009h for example has 
estimated ejecta masses for 5 type II SNe (1987A, 
1999em, 2003Z, 2005cs and 2004et) which would imply 
initial stellar masses a factor of ~2 higher than those 
estimated from the progenitor luminosity a nd stellar 
evolutionary tracks. lUtrobin fc Chugaii ( 20091) speculate 
that the hydrodynamics may be effected by asymme¬ 
try in the explosion and the lack of treatment of the 
multi-dimensional effects of Rayleigh-Taylor mixing be- 
twee n the h elium core a nd the hydrogen envelope. How¬ 
ever Spiro et al. ( 2014h used a different code (that of 
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Figure 7. From lO’Connor &: Ol~t] (1201 il l. The upper plot shows 
the compactness parameter £ 2.55 of stellar cores as a function of 
ZAMS mass for different stellar evolutionary models. The lower 
plot shows the time to black hole formation, assuming no ex¬ 
plosion has occurred. Reproduced from Figure 9 in “Black Hole 
Formation in Failing Core-Collapse Supernovae, by OConnor & 
Ott 2011, ApJ, 730, 70. 


Pumo & Zampieri 201 1) to find better agreement be¬ 


tween the direct progenitor masses and hydrodynamic 
masses estimates. A major goal for the SN community is 
to reconcile these independent and complimentary esti¬ 
mates as they fundamentally link the progenitor, explo¬ 
sion mechanism and SN observables. One would hope 
uniform agreement could be found that each method 
can inform the other of potential systematic errors. 


3.4 Black hole formation and failed 
supernovae 

On the basis of available data, the above two expla¬ 
nations are unlikely. Therefore the reason for the lack 
of high mass progenitors may be simply that they do 
not exist. In other words the SN population that we 
observe does not, on the whole, come from stars with 
initial masses greater than Mzams ~18M 0 . That these 
high mass stars will evolve to carbon-oxygen cores and 
subsequently to an Fe core composition and collapse 
is inevitable. However if black holes are formed, either 
directly or by fall back, and the neutrino energy depo¬ 
sition is not sufficient to drive the shock through the 
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steeper density profiles of these progenitors a failed SN 
is possible. It is plausible to postulate that the reason 
for the missing high mass progenitors is that they pro¬ 
duce failed SNe that have, so far gone undetected (or 
unrecognised) in all of our efforts to find and quan¬ 
tify transients in the local Universe. From a theoret¬ 
ical standpoint this is not surprising at all, the issue 
of explodability of the most massive stars has been a 
question since the 1980’s and 1990s. iFrver 


( 1999i) 


ex¬ 
plored black hole formation through direct collapse and 
fall back and already suggested then that stars above 
2OM 0 would produce black holes, possibly without a 
SN explosion. 

The density profiles in stellar cores is a strong func¬ 
tion of stellar mass, and recent work has revisited 
the idea that very massive stars may not produce 
detectab le SNe. As discuss ed in the introduction in 
( Sukhbold fe Wooslevl 2014 ) the density gradient in 
more massive stars is relatively shallow, which will likely 
lead to higher accretion rates onto the compact object 
created at the point of collapse. It would also lead to 
higher “ram pressure” which would need to be over¬ 
come by the shock (through the neutrino energy de¬ 
position mechanism) to launch a successful explosion 
through the star. The difficulty in getting successful 
simulated explosions has reinvigorated the study of the 
explodability of stars and in part i cular, three indepen¬ 
dent anal yses by O^Connor_& Ottl (120111 ). lu&iiano ct al 


( 2012 h ISukhbold fe Wooslevl (12014 ) and have used the 
compactness parameter (defined in O’Connor & Ott) 
which quantifies the core radius within which a certain 
amount of mass is confined in the final progenitor model 
£ = M/R(M) where M is in solar masses and R(M) is 
in units of 1000 km. Defining £ 2.5 for a value of 2.5M 0 
has been proposed. All three of these studies suggest 
that the compactness parameter is not a monotonically 
increasing and simple function of initial stellar mass. It 
will depend critically on mass-loss rate, internal mixing 
mechanisms in the stellar model not to mention metal- 
licity and rotation (and very likely mass-transfer in bi- 
narie s will have similar effect s to mass-loss). Figure[7]is 
from lO’Connor fe Ottl ( 201ll) showing the compactness 
C 2.5 at bounce, as a function of Mzams- The three stud¬ 
ies all suggest that at critical values of £ 2.5 > 0.2 — 0.3, 
successfid explosions are difficult to achieve. The figure 
illustrates that the compactness is dependent on the 
mass-loss recipes used in the stellar models, and that 
below Mzams < 2OM 0 , the stars are relatively easy to 
explode. However above that threshold, it becomes sig¬ 
nificantly more difficult to achieve explosions. There are 
what have been termed “islands of explodability” such 
that at some masses, the combination of mass-loss and 
final core structure may make it easier or hard to ex¬ 
plode. What is constant in all these three studies is 
that above Mzams > 2OM 0 , one would expect to have 
failed explosions and that there will be a tendency for 
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Figure 8. JJpper: The visually striking illustration of the disappearance of the red supergiant progenitor of SN2008bk, from 
iMattila et al.1 (|2010h . Panel A shows the VLT colour image of the progenitor (marked). Panel B shows the VLT NACO image of 
SN2008bk and the surrounding population at high resolution. Panel C shows an NTT colour image at approxim ately 940 days after ex - 
plosion illustrating the disappearance of the red source. The quantitative mass estimates of the pro genitor are i nlMaund et al.l <|2Q14ah . 
Middle and Lower: A spectrum of the faint blue source seen at 5 47 days after exp l osion from [Maguire et ah 1 te012P, showing the 
nebular phase emission lines from the SN ejecta, with a model from Jcrkstrand et al. ( 20^). The m odel is from an exploded 12Mq 
star, which is in good agreement with the updated progenitor mass estimate in lMaund et all |j2014ai l and in the summary table here 
(see Tablc[T] Reproduced from Figure 13, in “Constraining the physical properties of Type II-Plateau supemovae using nebular phase 
spectra”, K. Maguire et al, 2012, Monthly Notices of the Royal Astronomical Society, 420, 3451 


stars not to produce visible SNe. It may be that there 
is not a one to one relation and there is no simple, 
single numerical mass threshold at which neutron star 
or black hole formation occurs. Overall though, model 
stars above this threshold have a tendency to produce 
failed SNe at the point of core-collapse whereas stars be¬ 
low it generally proceed to explosion. In the same spirit, 


Horiuchi et al. i ( 2014t ) have carried out an extensive se¬ 


ries of simulations with more than 300 models and ar¬ 
gue that a critical value of £ 2.5 ~ 0.2 is the likely divide 
between successful and failed explosions. This would re¬ 
sult in stars in the initial mass range ~ 16 — 30Mq not 
producing canonical SNe. 
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The observational constraints presented in this review 
can reasonably be explained, and in a quantitative man¬ 
ner, with this model framework. The search for progeni¬ 
tors has been more difficult than was first thought, with 
the dearth of high mass and high luminosity progenitors 
an inescapable fact. Arguments can be made concern¬ 
ing dust, bias, and evolution to very faint WR stars but 
these have been addressed above. The reasonable agree¬ 
ment between the explosion theories based on the core 
compactness and the progenitor constraints mean they 
should b e considered to be phy sically linked. As first 
noted by iKochanek et al. ( 2008 1 the prediction of this 
scenario is that many massive stars (a fraction as high 
as / = 20 — 30%) end their lives by disappearing with¬ 
out a visible supernovae. As the dynamical timescale of 
a star of red supergiant dimensions is of order a year, 
these failed SNe may giv e rise to faint, long li ved re d 
transient s as p r opose d by Lovegrove fe Wooslev ( 2013ll . 
However Pirol ( 20131) proposed that a short timescale 
optical transient (lasting 3-10 days) with a faint lumi¬ 
nosity of ~ 10 40 — 10 41 ergs _1 could be produced due 
to the breakout of a shock induced in the stellar en¬ 


velope from neutrino losses in the core. Gerke et al 


( 20141) have just presented the first results from their 
novel and ongoing survey for failed stars and this has 
turned up one plausible candidate that requires further 
investigation. One could reasonably argue that black 
hole forming cores and failed explosions are not contro¬ 
versial ideas, rather they are solid theoretical predic¬ 
tions that hnve been__proposedjor more than 15 years 
(e.g. Fryer 19991) . Kochanek ( 2014 ) finds that the ob¬ 
served black hole mass function in the Galaxy can be 
modelled by assuming that black hole formation (and 
hence successful supernova explosion) is determined by 
the core compactness pa rameter . All o f th is adds im¬ 
petus to the searches of iGerke et al.1 ( 20141 ) and sur¬ 
veys for transients in the local Universe to identify the 
signatures t hat have been sugge sted by, for example, 
Lovegrove & Wooslev ( 20131) an d [Pirol (f2013i) . 

A final note on this topic of the explosion mechanism 
is the potential of recent work to compare the ejecta 
velocities in succes sful SNe t o the progenitor mass es¬ 
timates by Poznanski ( 120131 ). There appears to be a 
relation between the observed ejecta velocities, plateau 
duration and progenitor mass which could be explained 
if the energy deposition (the kinetic energy observed 
in the ejected envelope) is proportional to the initial 
mass of the progenitor cubed (E oc M 3 ). If this analy¬ 
sis holds then it could be telling us something about the 
efficiency of the deposition process in a neutrino revived 
explosion. As noted bv IPoznanskil ( 20131 ). the paucity of 
objects with low energy and large mass could support 
the idea of a population of failed SNe. This interesting 
correlation is worth pursuing further both from the ob¬ 
servational gathering of larger samples and theoretical 
investigation of the explosion mechanisms. 
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Figure 9. Upper Panel: This shows the observed line lumi¬ 
nosities for [Oi], Mgi] and Nai for SN2012aw compared t o the 
nebular models for 12 - 25Mq stars (Ijerkstrand et al]|2014l) . The 
line strengths support the mass estimate from direct detection 
of the progenitor star of around Mzams = 15M© and are much 
weaker than is predicted expect for exploding high mass pro¬ 
genitor models of 20-25Mq. Reproduced from Figure 4, in “ The 
nebular spectra of SN 2012aw and constraints on stellar nu¬ 
cleosynthesis from oxygen emission lines”, A. Jerkstrand et al, 
2014, Monthly Notices of the Royal Astronomical Society, 439, 
3649. Lower Panel: A compilation of the [Oi] line luminosi¬ 
ties compared to the predictions of exploding models of 12, 15, 
19 and 25M0. This again illustrates that type II SNe (predom¬ 
inantly II-P) do not have the observational signatures expected 
from higher mass 19-25M0progenitors. Reproduced from Figure 
9, in “Supersolar Ni/Fe production in the Type IIP SN 2012ec’\ 
Ijerkstrand et al. (l2015bl) Monthly Notices of the Royal Astro¬ 
nomical Society, 2015, 448, 2482 


3.5 Nebular spectra and nucleosynthesis 

A potentially powerful method of testing stellar evolu¬ 
tionary modelling and explosive nucleosynthesis is by 
observing, and physically modelling SNe when they are 
in their nebular phase. As the ejecta expand and be¬ 
come optically thin, the inner regions of the progen¬ 
itor star become visible and the elements synthesised 
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within the He and CO cores during stellar evolution 
become visible. Nebular spectra allow the diagnosis of 
nucleosynthesis products. This becomes particularly in¬ 
teresting for SNe with detected progenitors, as the inde¬ 
pendent estimate of the main sequence mass from neb¬ 
ular modelling can test the prediction from the progen¬ 
itor luminosity. An evolutionary model should repro¬ 
duce the pre-explosion stellar luminosity and the oxy¬ 
gen mass (for example ) prod u ced during helium and 
carbon burning, (as in Kozma &; Franssonlll998l) . Fig- 
ure[8] visually illustrates the diagnostic power of com¬ 
bining these two methods. The progenitor star is de¬ 
tected, seen to disappear and the faint blue flux visi¬ 
ble at the star’s position can be studied spectroscop¬ 
ically to determine element masses synthesised in the 
progenitor (e.g. the [O i] doublet at 6300,6364 is a use¬ 
ful diagnostic of the neutral oxygen mass). The key 
to this of course is having a reliable and sophisti¬ 
cated radiative transfer model that can reproduce the 
observed spectra and estimate element masses in the 
ejecta from the absolute and relative strengths of the 
emission lines. On e of the leading c odes for this has 
been develop ed bv l.Terkstrandl ( 201 ll i: [jerkstrand et al 


(2011 


Iop< 

12 ). 


2012). There are others such as Dessart et al 


( 2010 ) and iMazzali et ah ( 201Clh which have been suc¬ 


cessfully applied. The Jerkstrand et al. code, devel¬ 
oped in Stockholm, calculates the physical conditions 
in the expanding nebula and includes non-thermal heat¬ 
ing, ionization, and excitation from the gamma-ray and 
positron energy deposition of 56 Ni/ 56 Co (and 57 Co and 
44 Ti at later phases), and computes thermal and statis¬ 
tical equilibrium with the latest atomic data. It treats 
multi- line radiative transfer with a Monte Carlo tech¬ 
nique and is being actively applied to t hose SNe wi th 
progenitor d etections and limits (e.g. Ijerkstrand et al.1 
2014 . 2015alh These models can successfully reproduce 
the observed optical, near-IR and (when available) mid- 
IR spectra of ty pe II SNe at late phases (see the results 
for S N2004et in Jerkstrand et alJl2012t I 
20121 ) 


IMaguire et al 


In particular, the application to SN2012aw allows a 
self-consistent comparison of the progenitor luminos¬ 
ity and the oxygen produced. In a progenitor model, 
the stellar luminosity is determined by the He core lu¬ 
minosity. The He core size and burning rate dictates 
the oxygen production. Thus a progenitor model which 
has been artificially exploded should be able to re¬ 
produce the oxygen line strengths and the final lumi- 
nosity consistent l y. The physical sophistication of the 
Jerkstrand et al.l ( 201 lh code now allows this to be 
done to high accuracy. For SN2012aw the nebular phase 
spectra are in good agreement with the ejecta from a 
KEPLER model of M /ams = 15 M 0 progenitor star 
( Jerkstrand et al]l2014 ) and the forbidden neutral oxy¬ 
gen lines can be used to constrain the mass of oxygen 
to be less than < 1 M 0 (see Figure!!)]). This provides a 


physically consistent match to the progenitor luminos¬ 
ity of a KEPLER model of the same mass (14M 0 as in 
Tableland illustrated in FigureHJ). 

The [Ol] line strengths (and also the Mgi and Nal 
fluxes) are str ongly dependent on prog enitor mass. The 
conclusion of Ijerkstrand et ahl ( 2014 ) is that to date, 
there is no convincing example of a type II-P SN which 
displays the nucleosynthetic products expected from 
a Mzams > 20 M 0 progenit o r in it s nebular spectra 
(see Fig J9]) . Ijerkstrand et al.l ( 2015al) further modelled 
the nebular spectra of the lib SNe 1993J, 2008ax and 
201 ldh and constrained the initial mass of the progeni¬ 
tors to be Mzams = 12 — 16M 0 from the line strengths 
of the oxygen lines. Together with progenitor detections, 
constraints on the mass-loss rates, these nebular spectra 
and the radiative transfer analysis there is a consistent 
picture that there are no type II SNe in the local Uni¬ 
verse which come from high mass progenitors. While 
the nebular spectral analysis provides the strong pro¬ 
genitor constraints useful for this review it is becoming 
a powerful way to probe the enrichment of SNe and to 
allow quantitative checks on what role various SNe play 
in galactic chemical evolution models. 


3.6 Nucleosynthesis and chemical evolution 

If all stars above a certain main-sequence mass do 
form black holes then they are not likely to eject 
the elements synthesised during stellar evolution. The 
im portant implication s of t his were recently studied 
by iBrown &; Wooslevl ( 20131) who found that the so¬ 
lar abundance pattern (from oxygen to strontium, at 
A ~90) are fit well if there is no cut-off assumed. This 
calculation had all stars up to 12OM 0 produce suc¬ 
cessful super novae and produce nucle osynthetic yields 
according to Wooslev fe Hegerl ( 20071) . If the limiting 
main sequence mass that produces a black hole and no 
supernova (and therefore no enrichment) is reduced to 
Mbh = 25M 0 then the solar elemental abundance pat¬ 
tern is still comfortably reproduced. However there are 
unpalatable consequences if this mass limit is as low 
as Mbh = 18M 0 , for which there is some direct evi¬ 
dence as discussed earlier in this review. For this limit 
to be compatible with th e observed abundance pattern, 
Brown & Wooslev ( 20131) show that the mass-loss rates 
should be lower than currently assumed (to keep the 
carbon abundance under control), and the reaction rate 
for 22 Ne(a, n)Mg“ needs to be near the experimental 
maximum rate in order to reproduce the s—process iso¬ 
topes. In addition the historical supernova rate needs to 
be 3 times higher than that when no cut-off is assumed. 
The non-monotonic distribution of compactness dis¬ 
cussed in Section l3.4l which results in some higher mass 
stars exploding in the islands of explodability would 
mitigate these problems. However that would require 
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significant numbers of successful explosions above about 
22M 0 . One way out is that if supernovae from the most 
massive stars (M > 20Mq) are successful in enrichment 
but are more often enshrouded in dust and invisible to 
current surveys. The discussion in Section[3]3] disfavours 
this scenario. Searches for SNe in very high starfor- 
mation rate galaxies such as dusty luminous infra-red 
galaxies suggest that more than about 80% of expected 
core-collapse SNe are misse d in optical and near infra¬ 
red searches ( Mattila et al. 1 12013) . For normal galaxies, 
the missing fraction is more like ly to be 10-15% (again 
estimated in lMattila et al.f2012t ) and these missing SNe 
would need to be from the most massive stars to explain 
the deficit of high mass progenitors and ease the tension 
with the elemental production budget. Further work on 
searches for SNe specifically in the radio and mid-IR 
searches would be required to tie this down securely. 


4 CONCLUSION 

Drawing together the observational results by the com¬ 
munity over the period 1999-2013, together with the 
theoretical concepts on exploding stars, the following 
conclusions are proposed 


• There are 45 supernovae with either detected pro¬ 
genitors or upper limits in the volume limited sam¬ 
ple and either none (or one, if SN2009ip is in¬ 
cluded) has a progenitor above the luminosity limit 
of logL/L 0 ~ 5.1 (or an equivalent initial mass of 
M ~ 18Mq). A Salpeter IMF would suggest there 
should be 13 objects and the probability of finding 
0 or 1 is p — 3 x 10 -5 . 

• The various possible biases that have been sug¬ 
gested to explain this deficit were explored in this 
review including circumstellar dust, selection bias 
and luminosity analysis errors. It does not seem 
likely that they can explain the deficit. 

• One possible explanation is that most, or nearly 
all, stars above M /> 18M 0 evolve into WR stars 
that are hot and faint at the point of core-collapse 
when they produce Ibc SNe. However this ignores 
the growing evidence that the Ibc SN popula¬ 
tion are produced mostly by binary systems with 
masses in the range 8 — 2OM 0 . 

• Explosion modelling by several theoretical groups 
actually predict that stars above the mass limit of 
M ~ 18M 0 have a tendency to produce failed SNe 
and form black holes. Above this limit there are 
islands of explodability that can produce neutron 
stars and successful explosions. There is no mono¬ 
tonic relationship between the core “compactness” 
parameter (which will dictate the likelihood of a 
successful explosion) and initial mass. There is rea¬ 
sonable agreement between the observed limits for 


the luminosity of progenitors and the theory of ex¬ 
plosions. 

• The nebular spectra of core-collapse SNe which 
have identified progenitors is a new and powerful 
way to probe the progenitor mass range and nu¬ 
cleosynthesis in massive stars. The results to date 
support the idea that the progenitors discovered 
to date are in the 8 — 17M 0 range and illustrate 
that the supernovae from high mass stars which 
are thought to produce the bulk of cosmic oxygen 
have not yet been found. 
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